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We investigate exit times from domains of attraction for the mo- 
tion of a self-stabilized particle traveling in a geometric (potential 
type) landscape and perturbed by Brownian noise of small ampli- 
tude. Self-stabilization is the effect of including an ensemble-average 
attraction in addition to the usual state-dependent drift, where the 
particle is supposed to be suspended in a large population of identical 
ones. A Kramers' type law for the particle's exit from the potential's 
domains of attraction and a large deviations principle for the self- 
stabilizing diffusion are proved. It turns out that the exit law for the 
self-stabilizing diffusion coincides with the exit law of a potential dif- 
fusion without self-stabilization and a drift component perturbed by 
average attraction. We show that self-stabilization may substantially 
delay the exit from domains of attraction, and that the exit location 
may be completely different. 

1. Introduction. We examine the motion of a particle subject to three 
sources of forcing. First, it wanders in a landscape whose geometry is de- 
termined by a potential. Second, its trajectories are perturbed by Brownian 
noise of a small amplitude. The third source of forcing can be thought of as 
self-stabilization. Roughly, it characterizes the influence of a large popula- 
tion of identical particles subject to the same laws of motion. They act on 
the individual through an attractive potential averaged over the whole pop- 
ulation, which adds to the underlying potential drift. More formally, denote 
by Xf the random position of the particle at time t. It is governed by the 
d-dimensional SDE 

(1.1) dXf = V{Xf)dt- [ <5>{Xf -x)du £ t (x)dt + ^edW t . 



Received May 2006; revised October 2007. 

AMS 2000 subject classifications. Primary 60F10, 60H10; secondary 60K35, 37H10, 
82C22. 

Key words and phrases. Self-stabilization, diffusion, exit time, exit law, large devia- 
tions, interacting particle systems, domain of attraction, propagation of chaos. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2008, Vol. 18, No. 4, 1379-1423. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



S. HERRMANN, P. IMKELLER AND D. PEITHMANN 



In this equation, V denotes a vector field on M. d , which we think of as rep- 
resenting a potential gradient, the first source of forcing. Without the other 
two sources the motion of the particle would just amount to the dynamical 
system given by the ODE 

(1.2) ± = V(x). 

The small stochastic perturbation by Brownian noise W of intensity e ac- 
counts for the second source of forcing. It is responsible for random behavior 
of X s , and allows for transitions between otherwise energetically unreach- 
able domains of attraction. The integral term involving the process's own 
law u\ introduces a feature that we call self- stabilization. The distance be- 
tween the particle's instantaneous position Xf and a fixed point x in state 
space is weighed by means of a so-called interaction function <i> and inte- 
grated in x against the law of Xf itself. This effective additional drift can be 
seen as a measure for the average attractive force exerted on the particle by 
an independent copy of itself through the attraction potential <!>. In effect, 
this forcing makes the diffusion inertial and stabilizes its motion in certain 
regions of the state space. 

Equations of the type (1.1) are obtained as mesoscopic limits of microsys- 
tems of interacting particles, as the number of particles in an ensemble of 
identical ones tends to infinity, and subject to the same first two sources 
of forcing, that is, the force field V and the Brownian noise of intensity e. 
Suppose we are given an interaction function <I>, that is, for any two parti- 
cles located at x and y in state space the value <l>(x — y) expresses the force 
of mutual attraction. This attraction can, for instance, be thought of as 
being generated by electromagnetic effects. The dynamics of a particle sys- 
tem consisting of N such particles is described by the stochastic differential 
equation 

1 N 

dXt> N = V(Xl' N ) dt-±J2 $(Xl> N - X(> N ) dt + y/idWt, 

(1.3) 

X^ N = x l , l<i< N. 

Here the W % are independent Brownian motions. The self-stabilizing effect 
we are interested in originates in the global action of the system on the 
individual particle motion in the large particle limit N — > oo. Under suitable 
assumptions, in this limit the empirical measures YljLi ^ x j ' N can ^ e snown 
to converge to some law uf for each fixed time and noise intensity, and each 
individual particle's motion converges in probability to the solution of the 
diffusion equation 

(1.4) dXt = V(Xi)dt- f <S>{Xl-x)du e t {x)dt + ^dWl 
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The aim of this paper is to extend the well-known Kramers-Eyring law of 
exit from domains with noncritical boundaries by particles diffusing in po- 
tential landscapes with small Gaussian noise to systems (1.1) which include 
the described self-stabilization effect. In the potential gradient case without 
interaction, in which the individual particle's motion is interpreted by the 
solution trajectories Z e of the SDE 

(1.5) dZ £ t = -VU{Z e t ) dt + y/edW t , 

Kramers' law states that, in the small noise limit e — > 0, the asymptotic 
exit time of Z £ from a potential well of height H is of the order exp{^}. 
See the beginning of Section 4 for a precise formulation of this. We derive 
a similar statement for self-stabilizing diffusions. In particular we examine 
how self-stabilization adds inertia to the individual particle's motion, de- 
laying exit times from domains of attraction and altering exit locations. 
Mathematically, the natural framework for such an analysis is large devia- 
tions theory for diffusions. Our key ingredient for an understanding of the 
small noise asymptotics of the exit times proves to be a large deviations 
principle for self-stabilizing diffusions (1.1). In the potential gradient case, 
the rate function in the large deviations principle just minimizes the en- 
ergy needed to travel in the potential landscape. If the particle undergoes 
self-stabilization, energy has to be minimized in a landscape which addition- 
ally takes into account the potential of an attractive force that depends on 
the particle's distance from the corresponding deterministic path (1.2). Our 
main results (Theorems 3.4 and 4.2, 4.3) state that the large deviations and 
the exit behavior of X s are governed by this modified rate function. The 
techniques we employ to relate this time-inhomogeneous case to the classi- 
cal time-homogeneous one stipulate the assumption that the boundaries of 
the domains avoid critical points of the potential. 

Interacting particle systems such as (1.3) have been studied from various 
points of view. A survey about the general setting for interaction (under 
global Lipschitz and boundedness assumptions) may be found in [14]. There 
the convergence of the particle system to a self-stabilizing diffusion is de- 
scribed in the sense of a McKean-Vlasov limit, and asymptotic independence 
of the particles, known under the name propagation of chaos, as well as the 
link to Burgers' equation are established. Large deviations of the particle 
system from the McKean-Vlasov limit were investigated by Dawson and 
Gartner [4]. Further results about the convergence of the empirical distribu- 
tion of the particle system to the law of the self-stabilizing diffusion may be 
found in [3] or [10]. 

McKean studies a class of Markov processes that contains the solution of 
the limiting equation under global Lipschitz assumptions on the structure 
of the interaction [11]. A strictly local form of interaction was investigated 
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by Stroock and Varadhan in simplifying its functional description to a Dirac 
measure [13]. Oelschlager studies the particular case where interaction is 
represented by the derivative of the Dirac measure at zero [12]. Funaki ad- 
dresses existence and uniqueness for the martingale problem associated with 
self-stabilizing diffusions [9]. 

The behavior of self-stabilizing diffusions, in particular the convergence to 
invariant measures, was studied by various authors under different assump- 
tions on the structure of the interaction; see, for example, [1, 2, 15, 16]. 

The material in this paper is organized as follows. In Section 2 we discuss 
existence and uniqueness of strong solutions to (1.1). Strong solvability is 
nontrivial in our setting due to the self-stabilizing term, and is required for 
the subsequent investigation of large deviations. In Section 3 we derive and 
analyze the rate function modified by self-interaction, and this way obtain 
a large deviations principle for the diffusion (1.1). This proves to be the 
key ingredient for the analysis of exit times and a derivation of a version of 
Kramers' law for self-stabilizing diffusions in Section 4. We conclude with an 
illustration of our main results by discussing some examples which emphasize 
the influence of self-stabilization on exit time and exit location (Section 5). 

2. Existence and uniqueness of a strong solution. The derivation of a 
large deviations principle for the self-stabilizing diffusion (1.1) in the sub- 
sequent section involves pathwise comparisons between diffusions in order 
to apply the usual tools from large deviations theory, such as contraction 
principles and the concept of exponential equivalence. Their applicability 
relies on strong existence and uniqueness for (1.1), which is nontrivial in 
our situation since the solution process's own law appears in the equation. 
The interesting interaction term / &(X £ — x)duf(x) also adds a consider- 
able amount of complexity to the mathematical treatment. It depends on 
u\ = P o (Xf) ; thus classical existence and uniqueness results on SDE 
as well as the classical results on large deviations for diffusions (Freidlin- 
Wentzell theory) are not directly applicable. Consequently, the question of 
existence and uniqueness of solutions for (1.1) is an integral part in any 
discussion of the self-stabilizing diffusion's behavior, and will be addressed 
in this section. 

We follow Benachour et al. [1] to design a recursive procedure in order to 
prove the existence of the interaction drift b(t,x) = J $>(x — y)duf(y), the 
second drift component of (1.1). More precisely, we shall construct a locally 
Lipschitz drift term b(t, x) such that the classical SDE 

(2.1) dXf = V{X £ t ) dt - b{t, X £ t ) dt + VedWt, t > 0, 

admits a unique strong solution, which satisfies the additional condition 

(2.2) b(t,x)= [ $(x-y)d<(y)=E{$(x-X|)}. 
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In (2.1) W is a standard d-dimensional Brownian motion, and V :R d — >R d 
mimics the geometrical structure of a potential gradient. Of course, for (2.1) 
to make sense, the drift b has to be well defined, that is, the integral of (2.2) 
has to be finite, which depends upon certain moment conditions for X £ to 
be made precise later on. Apart from these moment conditions, existence 
and uniqueness for (1.1) will be understood in the sense that (2.1) and (2.2) 
hold with a unique b and a pathwise unique process X £ . 

For locally Lipschitz interaction functions of at most polynomial growth, 
Benachour et al. [1] have proved the existence of strong solutions in the 
one-dimensional situation, and in the absence of the vector field V. Since 
V forces the diffusion to spend even more time in bounded sets due to 
its dissipativity formulated below, it imposes no complications concerning 
questions of existence and uniqueness. Our arguments rely on a modification 
of their construction. 

Besides some Lipschitz type regularity conditions on the coefficients, we 
make assumptions concerning the geometry of V and 3> which render the 
system (3.1) dissipative in a suitable sense. All necessary conditions are 
summarized in the following assumption. 

Assumption 2.1. 

(i) The coefficients V and $ are locally Lipschitz, that is, for each R > 
there exists Kr > s.t. 

(2.3) \\V(x) - V(y)\\ + \Mx) - *(y)\\ < K R \\x - y\\ 

for x,y £ B R (0) = {z eR d :\\z\\ <R}. 

(ii) The interaction function is rotationally invariant, that is, there 
exists an increasing function (f> : [0, oo) — ► [0, oo) with (f)(0) = such that 

(2.4) <f>(x) = 1 ^(j)(\\x\\), xeR d . 

\\x\\ 

(iii) $ grows at most polynomially: there exist K > and r £ N such that 

(2.5) ||$(x)-$(y)|| < \\x-y\\(K+\\x\\ r + \\y\\ r ), x,y£R d . 

(iv) V is continuously differentiable. Let DV(x) denote the Jacobian of 
V. We assume that there exist Ky > and Rq > such that 

(2.6) (h,DV{x)h)<-K v 
for heR d s.t. \\h\\ = 1 and x 6 R d s.t. ||x|| > i? - 

The conditions that make our diffusion dissipative are (2.4) and (2.6). 
Equation (2.4) means that the interaction is essentially not more complicated 
than in the one-dimensional situation and has some important implications 
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for the geometry of the drift component E[$(x — Xf)] originating from self- 
interaction, namely that it points back to the origin. The same holds true for 
V due to (2.6). In the gradient case V = — VU, —DV is the Hessian of U, and 

(2.6) means that its eigenvalues are uniformly bounded from below (w.r.t. x) 
on neighborhoods of oo. Equation (2.5) is just a convenient way to combine 
polynomial growth and the local Lipschitz assumption in one condition. In 
the following two lemmas we summarize a few simple consequences of these 
assumptions. 

Lemma 2.2. There exist constants K,r],R± > such that the following 
hold true: 

(a) For all x,y£~R d 

(2.7) (x-y,V(x)-V(y))<K\\x-yf. 

(b) For x, y 6 M. d such that \\x — y\\ > FL\ 

(2.8) (x - y, V(x) - V{y)) < - V \\x - y\\ 2 . 

(c) For x£R d with \\x\\ > i?i 

(2.9) (x,V(x)} <-r]\\xf. 

Proof. Note first that, by continuity of DV, there exists K > such 
that 

{h,DV(x)h) < K 

holds for all x and all h of norm 1. Moreover, for x,y € M. d , x ^ y, we have 
V(x)-V(y) _ f 1 nT/ , , . x-y 



: / DV{y + t{x-y))^L.dt, 
Jo \\x-y\\ 



IF -2/11 Jo IF -2/1 

and therefore 

(2.10) /-£^iL /y ~ V l; y) )= f\h,DV(y + t\\x-y\\h)h)dt, 
\\\x-y\\ If -2/11 / •'0 

where h := iZ~y\\ ■ Since the integrand is bounded by K, this proves (a). 

For (b), observe that the proportion of the line connecting x and y that 
lies inside .Br o (0) is at most i u^m ■ Hence 

x-y V(x)-V(y)\ <K 2R Ky(l 2R ° 



\x — y\\ \\x — y\\ I \\x — y\\ \ \\x — y\ 

which yields (b). 
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Part (c) is shown in a similar way. Let x £ M. d with ||x|| > Rq, and set 
y := i?o]jfj[- Then the same argument shows the sharper bound 



/ x-y V(x)-V(y) \_/ x V(x)-V(y) \ 
XWx-yW' \\x-y\\ I \||x||' ||x||-J?o /' 



since the line connecting x and y does not intersect Br (0). Hence 



and 77=^. □ 

Lemma 2.3. For all x,y,z G M d we have: 

(a) \Mx - y )\\<2K + (K + 2 r + 1 )(\\x\\ r + 1 + \\y\\ r+1 ). 

(b) \\$(x- z )-<S>(y- z )\\ < \\ x -y\\[ K + 2 r (\\x\\ r + \\y\\ r + 2\\z\n}. 

(c) \\<S>(x-y)-<S>(x-z)\\<K 1 \\y-z\\(l + \\x\\ r ){l + \\y\\ r + \\z\\ r ), where 
Ki =max(K,2 r+1 ). 

(d) For all x, y £ R d and n £ N 

(2.11) (x\\x\\ n -y\\y\\ n ,$(x-y))>0. 

Proof. By (2.5) and since <1?(0) = we have 
\\Hx-y)\\<\\x-y\\(K + \\x-y\n 

<K(\\x\\ + \\y\\) + 2 r+1 (\\x\\ r+1 + \\y\\ r+1 ) 
< K(2 + \\x\\ r+1 + ||y|r +1 ) + 2'"+ 1 (||x|r+ 1 + ||y|r +1 ) 
= 2^+(^ + 2 r+1 )(||x|r +1 + ||y|r +1 ), 
that is, (a) is proved. For (b), we use (2.5) again to see that 

\\Q(x-z)-*(y-z)\\ < \\x-y\\(K+\\x-z\\ r + \\y-z\\ r ) 



Property (c) follows from x) = —<&(x) by further exploiting (b) as fol- 
lows. We have 



\ Mx _ y) _ Hx _ z) \l < || y _ 4[K + 2 r+1 (|| a ;|r + ||y|r + \\z\\ r )], 



which obviously yields (c) . Finally, (d) follows from a simple calculation and 
(2.4). Obviously, (2.11) is equivalent to (x||x|| n — ?/||y|| n ,rE — y) > 0. But this 
is an immediate consequence of the Schwarz inequality. □ 




Rv) + \\x\\\\V{y)\\ 



<iix- y ||[^+2'-(iixir + i| y ir+2||zir)]. 
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Let us now return to the construction of a solution to (1.1), that is, a 
solution to the pair (2.1) and (2.2). The crucial property of these coupled 
equations is that the drift b depends on (the law of) X s and therefore also 
on V, e and the initial condition xq. This means that a solution of (2.1) and 
(2.2) consists of a pair (X, b), a continuous stochastic process X and a drift 
term b, that satisfies these two equations. 

Our construction of such a pair (X, b) shall focus on the existence of the 
interaction drift b. It will be constructed as a fixed point in an appropriate 
function space such that the corresponding solution of (2.1) fulfills (2.2). Let 
us first derive some properties of b that follow from (2.2). 

Lemma 2.4. Let T > 0, and let (Xt)o<t<T be a stochastic process such 
that sup 0<t<T E,[\\Xt\\ r+1 ] < oo. Then b(t,x) = E[&(x — Xt)} has the following 
properties: 

(a) b is locally Lipschitz w.r.t. x G M. d , and the Lipschitz constant is 
independent of t G [0, T] . 

(b) (x - y, b(t, x) - b(t, y)) > for all x,y£R d ,t£ [0, T] . 

(c) b grows polynomially of order r + 1. 

PROOF. Note first that y \— > &(x — y) grows polynomially of order r + 1 
by Lemma 2.3(a), so that b is well defined. Moreover, we have 

||6(t, a0|| < E[\\<S>(x - X t )\\] <2K + {K + 2 r+1 )(||x|r +1 + E[\\X t \\ r+1 ]), 

which proves (c). For (a) observe that, by Lemma 2.3(b), we have for z £ 
x,y£B R (0), 

||$( x _ z ) _ $( y -z)\\< \\x - y\\[K + 2 r+1 (R r + \\z\\ r )]. 

Hence 

\\b(t, x) - b(t, y)\\ < E[\\<S>(x - X t ) - <% - X t )\\] 

<\\x-y\\[K + 2 r +\R r +E[\\X t f])] 

for x,y G Br(0). Since sup 0<t<T E[||X t || r+1 ] < oo, this implies (a). 
In order to prove (b), fix t G [0, T], and let n = Po Xf 1 . Then 

(x -y,b(t,x) -b(t,y)} 

\ x ~ y> TT — ~tMW x ~ U W) ~ TT — "ii^dly ~ n ID )fJ>(du). 
\ \\x — u\\ \\y — u\\ / 

The integrand is nonnegative. Indeed, it equals 

/ x — u 

\\x - u\\4>(\\x - u\\) + -u\\) - (y-u,- -<j>(\\x - u\\) 

\ \\x — u\\ 



d 
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- (x-U, Tp -rr 4>(\\y ~ U\ 

\ \\y- u \\ 



> \\x - u\\(f)(\\x - u\\) + \\y - u\\(f>{\\y - u\\) - \\y - u\\<p{\\x - u\\) 

- \\ x - u \\4>{\\y ~ u \\) 
= (\\x-u\\-\\y-u\\)(il>(\\x-u\\)-<f>(\\y-u\\)) 3 
which is nonnegative since (j) is increasing, so (b) is established. □ 

In the light of the preceding lemma it is reasonable to define a space of 
functions that satisfy the above stated conditions, and to look for a candidate 
for the drift function in this space. Let T > 0, and for a continuous function 
b:[0,T] xR d ^R d define 



(2.12) ||6|| r := sup sup 

te[o,T] x m d 1 + 



-V 



where q S N is a fixed constant such that 2q > r, the order of the polynomial 
growth of the interaction function <E>. Furthermore, let 

A T :={b:[0,T] x R d M. d \ \\b\\ T < oo,x \-+ b(t,x) 

(2.13) 

is locally Lipschitz, uniformly w.r.t. t\. 

Lemma 2.4 shows that, besides being an element of At, the drift of (2.1) 
must satisfy the dissipativity condition 

(2.14) (x-y,b(t,x)-b{t,y))>0, x,yeM d . 
Therefore, we define 

(2.15) A T :={beA T :b satisfies (2.14)}. 

It is obvious that || • ||t is indeed a norm on the vector space At- The 
subset At will be the object of interest for our construction of the interaction 
drift in what follows, that is, we shall construct the interaction drift as an 
element of At for a proper choice of the time horizon T. 

Once we have constructed the drift, the diffusion X will simply be given 
as the unique strong solution of (2.1) due to the following rather classical 
result about strong solvability of SDEs. It ensures the existence of a unique 
strong solution to (2.1) for a given drift b and is a consequence of Theorem 
10.2.2 in [13], since pathwise uniqueness, nonexplosion and weak solvability 
imply strong solvability. 

Proposition 2.5. Let (3 :R + xR d ^R d , (t,x) /3(t,x), be locally Lip- 
schitz, uniformly w.r.t. t € [0,T] for each T > 0, and assume that 

sup ||/?(t,0)|| < oo 

0<t<T 
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for all T > 0. Moreover, suppose that there exists ro > such that 
(x,P(t,x))<0 for \\x\\ > ro- 

Then the SDE 

dX t = l3(t,X t )dt + ^edW t 

admits a unique strong and nonexploding solution for any random initial 
condition Xq. 

It is easily seen that the drift P(t,x) = V(x) — b(t,x) does indeed satisfy 
the assumptions of Proposition 2.5 for any b £ At- This is an immediate 
consequence of (2.9) and (2.14). 

According to Lemma 2.4, the geometric assumptions formulated in As- 
sumption 2.1 imply that the drift term b of the self-stabilizing diffusion (1.1) 
is an element of At, provided the moment condition stated there is satis- 
fied. This moment condition is crucial for our construction of the drift, which 
motivates the following definition. 

Definition 2.6. Let T > 0. By a solution of (1.1) on the time interval 
[0,T] we mean a stochastic process (Xf)o<t<T that satisfies (2.1) and (2.2) 
on [0,T] and 

(2.16) sup E[||X t || 29 ] <oo. 

0<t<T 

A solution of (1.1) on [0,oo) is by definition a solution on [0,T] for each 
T>0. 

To construct a solution of (1.1) on [0,oo), we proceed in two steps. In the 
first and technically most demanding step, we construct a drift on a small 
time interval [0,T]. We shall define an operator T such that (2.2) translates 
into a fixed point property for this operator. To ensure the existence of a 
fixed point, one needs contraction properties of T which shall turn out to 
depend on the time horizon T. This way we obtain a drift defined on [0,T] 
such that the associated solution X exists up to time T. In a second step, 
we show that this solution's moments are uniformly bounded w.r.t. time, 
which guarantees nonexplosion and allows us to extend X to the whole time 
axis. 

To carry out this program, we start by comparing diffusions with different 
drift terms. 

Lemma 2.7. For b l ,b 2 £ At consider the associated diffusions 
dY t = V(Y t ) dt - b 1 ^, Y t ) dt + ^edW t 
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and 

dZ t = V{Z t ) dt - b 2 (t, Z t ) dt + y/edW t , 
and assume Yq = Zq. Then for t<T 

||F t -Z t ||<^ T ||6 1 -6 2 || T f\l + \\Z s \\ 2q )ds. 

Jo 

Proof. Since Y — Z is governed by a (pathwise) ODE, we have 
\\Y t - Z t \\ = " Z z s , V{Y S ) - V{Z S )^ ds 

^b l {s,Y s )-b l (s,Z s )Jds 

+ j^ ^ s i z z s sr b\s,z s )-b\s,z s )y s . 

The second integral in this decomposition is positive by definition of At, so 
it can be neglected. Furthermore, the first integral is bounded by K Jq \\Y s — 
Z s \\ds due to the dissipativity condition (2.7) on V. The last integral is 
bounded by 

f \\b 2 {s,Z s )-b 1 {s,Z s )\\ds<\\b 1 -b 2 \\ T f\l + \\Z s \\ 2q )ds. 
Jo Jo 

Combining these estimates yields 

\\Y t -Z t \\< K I \\Y S - Z s \\ ds + \\b l - b 2 \\ T f\l + \\Z s \\ 2q ) ds. 
Jo Jo 

Now an application of GronwalFs lemma completes the proof. □ 

The liberty of choice for the drift terms in Lemma 2.7 allows us to get 
bounds on Y and its moments by making a particular one for Z. We consider 
the special case of a linear drift term b(t,x) = Ax. 

Lemma 2.8. Let \ >K, and let Z be the solution of 

dZ t = V(Z t ) dt - XZtdt + VedW t . 

Furthermore, assume that E(||Zo|| 2m ) < oo for some m G N, m > 1. 
Then for allt>0 

E[||Z,|| 2m ]<2mt||y(0)||^- 1 exp{ £m((i+ fi 2 ^" 2)t } if Z = a.s., 
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and 

HI ^1 J - Ul UN j F| (E[||^o|| 2m l) 1/r " J 

+ 2mi||F(0)||^ m - 1 exp 

otherwise. 



em(d + 2m - 2)t 

w 



Proof. By Ito's formula we have for n > 2 

Z t \\ n = \\Z \\ n + M? + n [ \\Z s \\ n ~ 2 (Z s ,V(Z s )) - \\\Z s \\ n ds 

Jo 



(2.17) 



+ -(dn + n 2 -2n) / ||Z fl 
2 Jo 



where M n is the local martingale M f n = n^fe J^(Z s \\Z s \\ n - 2 ,dW s ). 
Since (x, V(x)) < —r]\\x\\ 2 for ||x|| > R\ according to (2.9), the first integrand 
of (2.17) is negative if \\Z S \\ > R\- If \\Z S \\ < we use the global estimate 
{x,V(x)) < K\\x\\ 2 + 11^(0)1111x11, which follows from (2.7). We deduce that, 
since A > K, 



- 2 (z s , v{z s )) - \\\z s \\ n <{k- \)\\z s \\ n + ||^(o)|| ||z.r 

n— 1 



<\\vm\R n i 

Thus, 

\\Z t \\ n < \\Z \\ n + M? + n\\V(0)\\tRi~ 1 + Udn + n 2 - 2n) f \\Z, 

2 Jo 

Using a localization argument and monotone convergence yields 

EOlZtlH^EOI^bH+nllFCO)!!^- 1 
(2.18) t 

+ -(dn + n 2 - 2n) f E[\\Z s \\ n ~ 2 } ds . 
2 Jo 

We claim that this implies 



2m l <X^w\\\ zJI 2 ( m -i)l^lL 



E[\\z t \\ 2m ]<j2n\\M 2{m ' j) } 



(2.19) 



i=o 



J! 



for all m £ N, m > 1, where ce m = em(d + 2m — 2). Indeed, for m = 1 this 
is evidently true by (2.18). The general case follows by induction. Assume 
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(2.19) holds true for m- 1. Then by (2.18) 

E[||Z t || 2m ] <E[||Z || 2m ]+2m||F(0)||t J R? m - 1 

rt — - ",(a m _is) J_1 



+ a m / ^E[||Zo|| 2 ( m -^]- 
Jo ~1 



+2(m - 1) gffl flr -f : is 

< E[||Zo|| 2 '"] + 2m||V(0)||«? m - 1 

+ ^a m E[||Z || 2 ( m -^]^^ 
i=i ] ' 

m oP~ 2 P 
+ 2m\\V(0)\\Ri m - 1 Y, a m^\y- 
3=2 Rl ?■ 
m j ,j 

< 2m||T/(0)pi? 2m - 1 +^E[||Z || 2(m - i) ]^p- 

m j—i ,j 

+ 2m\\V(0)\\Rl m+1 T,%T 

3=1 K l J- 

m j ,j m j—l-ij 

= ^E[||Z || 2 ^)]^i- + 2m\\V(0)\\Ri m+1 E ^JT' 

3=0 j=l Rl J- 

and so (2.19) is established. Since E[||Z || 2(m ^] < (E[||Z || 2m ]) 1 - J ' /m for 
j < m, we may exploit (2.19) further to conclude that 

m j ,j 

E[||Z t || 2m ]<E[||Z || 2m ]E 



f^j\(E[\\Zo\\ 2m ]V /m 



+ 2mt\\V(0)\\Ri m - 1 J2 
<E[||Z || 2m ]exp 



3=1 K l J ! 



E[\\z \\ 2m ]y/r, 



2mt||y(0)||i? 2m - 1 exp 



Rl { 



which is the announced bound if we identify the first term as zero in case Zq = 0. 
□ 
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Let us define the mapping T on At that will be a contraction under 
suitable conditions. For fr£ At, denote by X^ the solution of 

(2.20) dX t = V(X t ) dt - b(t, X t ) dt + ^edW t , 

and let Tb(t,x) :=E[$(x — xj®)]. By combining the two previous lemmas, 
we obtain the following a priori bound on the moments of X^ b \ 

Lemma 2.9. If the initial datum of (2.20) satisfies E[||xJ b) \\ 2qn ] < oo for 
some n £ N, then for each T > there exists k = k(n, T) > such that for 
all be A T 

sup E[||X t (6) |n < k{l + T n e nKT {\\b\\ T + K n )). 

0<t<T 

Proof. Let b l (t,x) := b(t,x) and b 2 (t,x) = Kx, and denote by Y, Z the 
diffusions associated with b , b 2 . By Lemma 2.7 we have for t £ [0,T] 

E[||y t n < 2"(E[iiz t ir i ] +E[||y t - Ztin) 

< 2™E[||Zi|n + 2 n e nKT t n ~ 1 \\b 1 - 6 2 ||^e[ f\l + ||Zj 2<? ) n (is 

Uo 

<2 n (l + E[||Zt|| 2gn ]) 

+ 2 n e nXT t n (||6 1 || T + ||6 2 || T ) n sup E[(l + ||Z s || 2<? ) n ] 

0<s<T 

2 <?™i \ n i ^"^ni^iia _l ii a 2 nil 1 



<8 n 1+ sup E[||Z s ||^ n ] )(l + t n e nKI (\\b l \\ T + \\b Z \\ T )). 

\ 0<s<T J 

Due to the assumption E[||X^ 6) \\ 2qn ] < oo, the constant k(n, T) = 8 n (l + 
sup 0<s<T E[||Z s || 2ljn ]) is finite by Lemma 2.8. Furthermore, we have ||6 2 ||t < 
K, that is, the lemma is proved. □ 

Now we are in a position to establish the local Lipschitz continuity of the 
operator V. The explicit expression for the Lipschitz constant shows that T 
will be a contraction on a sufficiently small time interval. 

Lemma 2.10. Let b x ,b 2 £ At, and denote by Y,Z the corresponding 
diffusions as in Lemma 2.7. For i £ N let rrii(T) = sup <£<rE[||Y t ||*] and 
n 4 (T)=sup <,< r E[||^|H. 

There exists a constant k = k(m4 q (T),n4 q {T)) such that 



ir&i - r& 2 ||T < kVTe^Wb 1 - b 



2 1 
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Proof. From Lemma 2.3(c) and the Cauchy-Schwarz inequality follows 
that 

\\Tb\t,x)-Tb 2 (t,x)\\ 

<E[\mx-Y t )-<t>(x-Z t )\\] 

<k 1 (i+ iixir)E[i|y* - z t \\(i + \\Y t r + ii^id] 



<K 1 (1 + \\x\\ r )y/E[\\Y t - Z t m(l + \\Y t \\r + \\Z t \\n 2 }, 

where K\ = max(i^, 2 r+1 ). By Lemma 2.7, since (1 + x) 2 < 2(1 + x 2 ), we 
have 

E[||y - Z t \\ 2 } < e 2KT \\b x - b 2 \\ 2 T E^j\l + ||^ s || 2 ^)^ 2 " 

< e 2^T ||6 l_ 62|| ^ [ T m + \\Z S \\ 2 1) 2 } ds 

Jo 

<2Te 2KT \\b 1 -b 2 \\ 2 T (l+ sup n\\Zs\\ 4g }). 

\ 0<s<T / 

Moreover, using the inequality (a + b) 2 < 2(a 2 + b 2 ), we deduce that 
E[(l + ||y t f + \\Z t \\ r ) 2 } < 2(1 + 2E[||y \\ 2r + \\Z t \\ 2r ]) 

<io(i + E[||y t || 4 « + ||z t || 4 «]), 

where we exploited that 2q > r implies E[\\Y t \\ 2r ] < 1+E[||y || 4 «], and likewise 
for the moment of Zt- By combining all these estimates, we find that 

\\Yb\t,x)-Tb 2 (t,x)\\ 
1 + ||x|| 2 9 

l + iizir 



<2K 1 V5Te KT \\b 1 -b 2 \\ T 



1+ ||x|| 2( 7 



1/2 

x (i+ sup e[\\z s \\^}) (i + EHiyii^ + 11^11^]) 1 / 2 . 



0<s<T 

Hence, if we set k := 4i^i\/5{(l + n^ q {T)){l + rrn q {T) +n^ q (T))} 1 ^ 2 , we may 
conclude that 

\\Tb 1 -rb 2 \\ T <kVTe KT \\b 1 -b 2 \\ T , 
that is, k is the desired constant. □ 

The next proposition shows that the restriction of T to a suitable subset of 
the function space At is a contractive mapping, which allows us to construct 
a solution on a small time interval. 
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Proposition 2.11. Forv > let A T = {b e A T : ||6||t < Assume i/iai 
the initial condition Xq satisfies E[||Ao|| 29n ] < oo for some n > 4(7. There 
exists vq > suc/i £/ia£ for any v > vq there exists T = T(v) > such that 
the following hold true: 

(a) r(Ay) C A T , and £/ie Lipschitz constant ofT\A T is less i/ian i. 

(b) There exists a strong solution to (2.1), (2.2) on [0, T] which satisfies 

sup E[||X t (b) || n ] < yfc(l + T n e n ^ T (^ n + i^ n )), 

0<t<T 

where k = k(n,T) is the constant introduced in Lemma 2.9. 

Proof. Let b e A T , and let X = X® and m^T) = sup < t < T £[11^11*] 
for i £ N. By Lemma 2.9 the condition E[||Xo|| 2ljn ] < oo implies rai(T) < oo 
for T > and i <n. Moreover, Lemma 2.3 shows that 

||r&(t,a;)|| < 2K + (K + 2 r+l )(\\x\\ r+1 +E[\\X\\ r+l ]) 

<K(l + \\x\\ r+1 )(l + E[\\X t \\ r+1 ]), 

where K = 2K + 2 r+1 . Consequently, by definition of || • \\t, 

(2.21) ||r6|| T <2J?(l + m r+ i(T)), t<T. 
By Lemma 2.9 there exists k = k(r + 1, T) > such that 

(2.22) m r+1 (T) < k(l + r r+1 e (r+1) ^ r (||6||5, +1 + K r+1 )). 

This inequality, together with (2.21), is the key for finding a suitable subset 
of At on which T is contractive. The r.h.s. of (2.22) converges to k as T — > 0, 
and this convergence is uniform w.r.t. 6 £ A^ for each > 0. The dependence 
of the limiting constant k on T imposes no problem here; just fix k = k(r + 
l,To) > for some To and use the fact that (2.22) is valid for all T < Tq, as 
the proof of Lemma 2.9 shows. 

Thus, we may fix i/q > 2K(1 + k) and deduce that for any v > uq we can 
find To = Tq(v) such that ||6||t < v implies ||r&||T < v for T < Tq. Moreover, 
by Lemma 2.4, Tb satisfies all the conditions as required for it to belong to 
At, that is, T maps A T into itself for all T <Tq. Additionally, the assump- 
tion n> 4q implies that m^T) is uniformly bounded for all b in A T , and 
Lemma 2.10 shows that, by eventually decreasing To, we can achieve that 
r is a contraction on A T with Lipschitz constant less than ^, that is, (a) is 
established. 

In order to prove (b) , the existence of a strong solution on the time interval 
[0, T] for some T < Tq, we iterate the drift through T. Let 6q G A t , and define 



bi + \:=Tbi for is Nq. 
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The contraction property of T yields ||6j+i — hllx <2 l ||&i — 6o||t for all i, 
and therefore 

oo 
n=0 

which entails that (bi) is a Cauchy sequence w.r.t. || • By definition of 
|| • \\t, (pi) converges pointwise to a continuous function b = b(t,x) with 
||6||t < 00. It remains to verify that the limit is again an element of At- In 
order to see that it is locally Lipschitz, let := X^ bi \ As in the proof of 
Lemma 2.4, we have for x,y G -Br(O) 

WTbtfrx) - rbi(t, y)\\ < E[||d>(* - X® ) - *(y - X«)||] 

<\\x-y\\[K + 2 r+1 (R r +E[\\X®\\ r })}. 

Since ||&i||T < v for all i, (2.22) yields 

sup sup E[\\xf\\ r ]<k(l + T r+1 e^ KT (v r+l + K r+l )). 

iGN 0<i<T 

Therefore, we may send i — > 00 to conclude that 6 is locally Lipschitz. b being 
the pointwise limit of the bi, it inherits the polynomial growth property and 
the dissipativity condition as stated in Lemma 2.4(b) and (c). (Notice that 
we may not invoke Lemma 2.4 at this stage.) 

It remains to show that the diffusion X = X^ associated to b has the 
desired properties. Note first that the existence of X is guaranteed by the 
classical result of Proposition 2.5. Since Tb = b, which means that 

b(t,x) =Fb(t,x) =E[$(x- Af°)] 

for t £ [0,T] and x G M. d , X is the diffusion with interaction drift b. The 
boundedness of its moments is again a consequence of Lemma 2.9. □ 

Let us recall the essentials of the construction carried out so far. We have 
shown the existence of a solution to (1.1) on a small time interval [0, T]. For 
the moments of order n to be finite, one needs integrability of order 2qn 
for the initial condition. Moreover, the parameter n needs to be larger than 
or equal to Aq in order for the fixed point argument of proposition 2.11 to 
work. Observe that the condition n> 4q appears first in this proposition, 
since this is the first time the process is coupled to its own drift, while 
in all previous statements the finiteness of moments is guaranteed by the 
comparison against the diffusion Z, which is governed by a linear drift term. 
In order to find a solution that exists for all times, we need to carefully 
extend the constructed pair (X,b) beyond the time horizon T. Although 
nonexplosion and finiteness of moments would be guaranteed for all T by 
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Proposition 2.5 and Lemma 2.9, we have to take care of the fact that the 
drift itself is defined only on the time interval [0,T]. With sufficiently strong 
integrability assumptions for Xq one could perform the same construction on 
the time intervals [T,2T], [2T, 3T] and so on, but one loses an integrability 
order 2q in each time step of length T. 

For that reason we need better control of the moments of X over the 
whole time axis, which is achieved by the following a posteriori estimate. 

Proposition 2.12. Let m G N, m> 4q 2 , such that E[||X || 2m ] < oo. 
For each n 6 {1, . . . , m} there exists a constant a = a{n) > such that the 
following holds true for all T > 0: if X solves (1.1) on [0, T], then 

sup E[\\X t \\ 2n ]<a(n). 
0<t<T 

Proof. Fix T > 0, and assume that X solves (1.1) on [0, T] (in the sense 
of Definition 2.6). Then b(t, x) := E[$(x - X t )} belongs to k T by Lemma 2.4, 
and m > Aq 2 implies sup <t<TlE[||Xj|| 4 ' ? ] < oo by Lemma 2.9. Let / n (t) = 
E[||Xt || 2n ]. We proceed in several steps. 

Step 1: Boundedness in L 2 . We know already (by Lemma 2.9) that 

sup fi{t) < oo. 

0<t<T 

The only point is to show that the bound may be chosen independent of T. 
By Ito's formula we have 

/i(t) = E[||X || 2 ]+rtd + 2 [ E[{X s ,V(X s ))]ds-2 [ E[(X s ,b{s,X s ))]ds. 

Jo Jo 

Let us estimate the last term that contains the interaction drift b. Note first 
that 2q> r implies r + 2 < 4g, so sup < t <r E[||X t || r+2 ] < oo as pointed out at 
the beginning of the proof, and the Cauchy-Schwarz inequality yields that 
M[(Xt,b(t,Xt)}] is finite for t £ [0, T] since b grows polynomially of order 
r + 1. By definition of 6, we may take an independent copy X of X, to write 

2E[(X S! b(s, X s ))] = 2E[(X S , <S>(X S - X s ))] 

= E[(X S , $(X S - X s ))} - E[{X S , $(X S - X,))] 

= E[{X s -X s ,<S>(X s -X s ))}>0 

where the last inequality is due to (2.4). In order to estimate the other 
integral, let R> R±. Using (2.9) and the local Lipschitz property of V, we 
see that 

E[{X s ,V(X s ))]<-r,E[\\X s \\ 2 l {llx4>R} ] 

+ E[(K\\X s \\ 2 + \\V(0)\\\\X s \\)l {llXs \\<R}] 
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< -??E[||X S || 2 ] + ( V + K)R 2 + \\V{Q)\\R 

= - r] f 1 (s) + R{\\V{Q)\\+R(r ] + K)). 

Obviously, f\ is differentiable, and summing up these bounds yields 

f[(t) <ed- 277/1 (t) + 2i?(||F(0)|| + R( V + K)). 

Thus, there exists 7 > such that {t G [0, T] : fi(t) > 7} C {t G [0, T] : f[(t) < 
0}, which implies fi(t) < /i(0) V7 for all t G [0, T]. This is the claimed bound, 
since 7 is independent of T. 

Step 2: Moment bound for the convolution. Let X be an independent 
copy of X, that is, a solution of (1.1) driven by a Brownian motion that 
is independent of W. In this step we shall prove that E[||A^ — A(|| 2n ] is 
uniformly bounded w.r.t. time. 

Let R > R lf and let r = inf{t > 0:\\X t - X t \\ > R}, g n (t) = E[\\X t - 
X t || 2 "l {t<r} ] and w n (t)=E[\\X tAT -X TAt \\ 2n ]. Then w n (t) = g n {t) + R 2n F(t > 

t). Furthermore, using the SDE (1.1) for both X and X, applying Ito's for- 
mula to the difference and taking expectations, we obtain for n > 1 

rtAT 



w n (t) = E[\\X - X \\ 2n ] + en(d + 2n- 2)E 



\Xf, — X s 



\2n-2 



ds 



+ 2nE 



2nE 



Mr 



Mr 



\X X — X s 



\X S — X K 



\2n-2 



\2n-2 



(X s -X s ,V(X s )-V(X s ))ds 



(X s -X s ,b(s,X s )-b(s,X s ))ds 



The last term is negative by Lemma 2.4, which yields together with (2.7), 
(2.8) and Holder's inequality 

w' n {t) < en(d + 2n- 2)E[\\X t - X t \\ 2n - 2 l {t<T} ] 

+ 2nE[\\X t - X t \\ 2n - 2 (X t - X t ,V(X t ) - V(X t ))l {t<r} ] 



< en(d + 2n - 2)g n _ 1 {t) 
+ 2n(K + r])E[\\Xt - 



X t \\ 2n l 



{\\Xt-Xt\\<Ri;r>t}\ 



2nr,E[\\X t -X t \\ 2n l {t<T} ] 



< en(d + 2n - 2)g n (t) 1 ~ 1/n + 2n(K + r])R\ n - 2nr)g n (t). 

As in the first step, there exists some constant 5 > such that {t G [0, T] : g n (t) > 
5} C {t G [0,T]:w' n (t) < 0}. Since w n — g n is nondecreasing this implies 
gn{t) < 5n(0) V <5 for all t G [0, T]. Moreover, <5 depends only on the con- 
stants appearing in the last inequality and is independent of the localization 
parameter R. Hence, by monotone convergence, we have 



E[||X* - X t \\ 2n ] < E[\\X - X \i zn ] V5, t G [0,T]. 



2rn 
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Step 3: Bound for the centered moments of X. In this step we shall prove 
that the moments of Yt := Xt — E[Xt\ are uniformly bounded. We proceed 
by induction. The second moments of X are uniformly bounded by the first 
step; so are those of Y . Assume the moments of order 2n are uniformly 
bounded by j n > 0. If n + 1 < m, we may invoke step 2, to find <5 n +i > 
such that E[\\X t - X t \\ 2n+2 ] < 5 n+1 for t £ [0,T]. Now we make the following 
observation. If £, £ are independent, real-valued copies of each other with 
E[f] = 0, then 

E[(f - i) 2n + 2 ] = 2E[£ 2n+2 } + ^ ( 2n + 2 ] (-l) k E[£ k }E[f n+2 - k ], 

k=2 ^ ' 

and therefore 

2n 

2E[£ 2n+2 ] < E[(f - |) 2n+2 ] + ( 2n ^ 2 ) |E[^]E[e 2n+2_fc ]| 

< e[(^ - £) 2 ™+ 2 ] + 2 2n+2 (l + E[£ 2n ]) 2 . 

Let us apply this to the components of Y, and denote them by Y 1 , . . . , Y d . 
We obtain for t G [0, T] 



2E[||yf|| 2n+2 ] < 2d n+1 E 



d 

\2n+2 



Lj=l 



< d n+1 J2E[(X° t - X|) 2n + 2 ] + 2 2n+2 (l + E[(y/) 2n ]) 2 

3=1 

< d n+2 (E[||X t - X t || 2n + 2 ] + 2 2n+2 (l + E[||Y t || 2n ]) 2 ) 

<^ +2 (,5 n+1 + 2 2 ™ +2 (l + 7n) 2 ), 

which is a uniform bound for the order 2{n + 1). 

Step 4: Bound for the moments of X. In the fourth and final step, we prove 
the announced uniform bound for the moments of X. It follows immediately 
from the inequality 

E[\\X t \\ 2n } < 2 2n (E[\\X t - E[X t ]\\ 2n ] + ||E[X t ]|| 2 ™). 

The last term satisfies ||ELYj]|| 2n < /i(t) n , which is uniformly bounded ac- 
cording to step 1, and the centered moments of order 2n are uniformly 
bounded by step 3 whenever n<m. □ 



The results concerning the existence of X s are summarized in the follow- 
ing theorem. 
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Theorem 2.13. Let q := [§ + 1], and Zei A~o be a random initial con- 
dition such that E[||Xo|| 8lj2 ] < oo. Then there exists a drift term b(t,x) = 
b £,x °(t,x) such that (2.1) admits a unique strong solution X £ that satisfies 
(2.2), and X £ is the unique strong solution of (1.1). Moreover, we have for 
allneN 

(2.23) supE[||Xf|| 2n ] <oo 

whenever E[||.Xg|| 2n ] < oo. In particular, if Xq is deterministic, then X £ is 
bounded in L P (F <g) Arg^i) for all p > 1. A is used as a symbol for Lebesgue 
measure throughout. 

Proof. In a first step, we prove uniqueness on a small time interval. 
Let K = 2K + 2 r+1 , and choose a(q) > according to Proposition 2.12. 
By Proposition 2.11 there exist v > 2K(2 + a(q)), T = T(u) > and b G 
A r such that Tb = b, that is, X = X^> is a strong solution of (1.1) on 
[0,T]. Assume Y is another solution of (1.1) on [0,T] starting at Xq such 
that m 2q {T) := sup <i< T E[||y 4 || 29 ] < oo, and let c(t,x) =E[$(x-F t )]. Then 
c G At by Lemma 2.4, and Tc = c. Moreover, it follows from (2.21) and 
Proposition 2.12 that 

\\c\\ T < 2K(2 + m 2q (T)) < 2K(2 + a(qf) < u, 

that is, c G A T . Hence c is the unique fixed point of T|A T . Thus c = b, and 
Proposition 2.5 yields X = Y. 

In the second step, we show the existence of a unique solution on [0, oo). Let 
U := sup|t > : (1.1) admits a unique strong solution X on [0, T], 

sup E[||A: t || 2 ' ? ] <oo 

0<t<T 

By the first step we know that U > 0. Assume U < oo. As in the first 
step, choose a(4q 2 ) > according to Proposition 2.12, and then fix z> > 
2if (2 + a{4q 2 )) and T = f(v) > that satisfy Proposition 2.11. Let < 5 < 
mm(U, T/2), and fix T £]U — 8, U\. There exists a unique strong solution X 
on [0,T], and E[||X T || 8,?2 ] < oo by Proposition 2.12. Now consider (1.1) on 
[T, oo) with initial datum Xt- As in the first step, we may find a unique 
strong solution on [T,T + T]. But this is a contradiction since T + T > U. 
Consequently, U = oo, and (2.23) holds by Proposition 2.12. □ 
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3. Large deviations. Let us now turn to the large deviations behavior of 
the diffusion X s given by the SDE (1.1), that is, 

(3.1) dX\ = V(Xf) dt - f -x)du £ t (x)dt + y/edW t ,t>0, 

Xq = Xq £ R d . 

The heuristics underlying large deviations theory is to identify a determin- 
istic path around which the diffusion is concentrated with overwhelming 
probability, so that the stochastic motion can be seen as a small random 
perturbation of this deterministic path. This means in particular that the 
law uf of Xf is close to some Dirac mass if e is small. We therefore pro- 
ceed in two steps toward the aim of proving a large deviations principle for 
X s . In a first step we "guess" the deterministic limit around which X s is 
concentrated for small e, and replace uf by its suspected limit, that is, we 
approximate the law of X s . This way we circumvent the difficulty of the 
dependence on the law of X s — the self-interaction term — and obtain a dif- 
fusion which is defined by means of a classical SDE. We then prove in the 
second step that this diffusion is exponentially equivalent to X s , that is, it 
has the same large deviations behavior. This involves pathwise comparisons. 

3.1. Small noise asymptotics of the interaction drift. The limiting be- 
havior of the diffusion X s can be guessed in the following way. As explained, 
the laws u\ should tend to a Dirac measure in the small noise limit, and since 
$(0) = the interaction term will vanish in the limiting equation. Therefore, 
the diffusion X £ is a small random perturbation of the deterministic motion 
i/j, given as the solution of the deterministic equation 

(3.2) ip t = V(i/} t ), i>o = x , 

and the large deviations principle will describe the asymptotic deviation of 
X s from this path. Much like in the case of gradient type systems, the 
dissipativity condition (2.9) guarantees nonexplosion of tp. Indeed, since 
^HV'tll 2 = 2(tj)t,ipt) = 2(A,V{i/j t )}, the derivative of ||^|| 2 is negative for 
large values of ||V>t|| by (2.9), so ip is bounded. In the sequel we shall write 
ipt{xo) if we want to stress the dependence on the initial condition. 

We have to control the diffusion's deviation from this deterministic limit 
on a finite time interval. An a priori estimate is provided by the following 
lemma, which gives an L 2 -bound for this deviation. For notational conve- 
nience, we suppress the e-dependence of the diffusion in the sequel, but keep 
in mind that all processes depend on e. 

Lemma 3.1. Let Z t := X t — tpt(x ). Then 

E\\Z t \\ 2 <etde 2Kt , 
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where K is the constant introduced in Lemma 2.2. In particular, Z — > as 
e — > in L P (P (g) Ap^l) f or all p>l and T > 0. 77ns convergence is locally 
uniform w.r.t. the initial condition xq. 

Proof. By Ito's formula we have 

\\Z t \\ 2 = 2^I f\z s ,dW s )-2 f t (Z s ,b £ ' X0 (s,Z s + ^ s {x )))d S 
Jo Jo 

+ 2 f\z s , V(Z S + tP s (x q )) - V(ip s (x ))) ds + etd. 
Jo 

By Theorem 2.13 X and thus Z is integrable of all orders. In particular, Z is 
square-integrable, so the stochastic integral in this equation is a martingale. 
Now consider the second term containing the interaction drift b £ ' x ° . Let v s = 
Po Z~ l denote the law of Z s . Since Z has finite moments of all orders, Lemma 
2.3 implies / J \\(z,$(z — y))\\v s (dy)v s (dz) < oo. Thus, by Assumption 2.1(h) 
about the interaction function $ and Fubini's theorem, 



2E{Z s ,b £ ' xo (s,Z s + Mxo))) = 2 I (z,E[*(z + M*o)-X')])M<b) 

{z,®(z - y))v s {dy)v s (dz) 



(z - y,$(z - y))u s {dy)v s (dz) > 0. 
Hence by the growth condition (2.7) for V 

E\\Z t \\ 2 < 2 f l E(Z s ,V(Z s + ^(xo)) - V(t/,.(xo))) ds + etd 



o 



<2K I E\\Z s \\ 2 ds + etd, 
o 



and Gronwall's lemma yields 

E\\Z t \\ 2 <etde 2Kt . 

This is the claimed bound. For the L p -convergence observe that this bound 
is independent of the initial condition xq. Moreover, the argument of Propo- 
sition 2.12 shows that sup{E(||X t || p ) : < t < T,x £ L,0 < e < e } < oo 
holds for compact sets L and eo > 0. This implies that Z is bounded in 
L p (P<g> Ar 0) Ti) as e — > 0, uniformly w.r.t. xo G L. Now the L p -convergence 
follows from the Vitali convergence theorem. □ 



Corollary 3.2. For any T > we have 

lim b £ ' x °(t,x) = <E>(x - ip t (xo)), 

s~^0 

uniformly w.r.t. t G [0, T] and w.r.t. x and xq on compact subsets ofM. d . 
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Proof. The growth condition on «3> and the Cauchy-Schwarz inequality 
yield 

\\b £ (t,x)-<f>(x-Mxo))\\ 2 

< E[\\X t - Mxo)\\(K+\\x - X t \\ r + \\x - V^(a=o)ir)] 2 

< E[\\x t - Mxo)\\ 2 m(K + ii* - x t \\ r + ii* - Vt(^o)ir) 2 ]. 

The first expectation on the r.h.s. of this inequality tends to zero by Lemma 
3.1. Since X is bounded in L 2r (F), uniformly w.r.t. xq on compact sets, the 
claimed convergence follows. □ 

In a next step we replace the diffusion's law in (3.1) by its limit, the Dirac 
measure in ip t (xo). Before doing so, let us introduce a slight generalization 
of X. 

Theorem 2.13 implies that X is a time-inhomogeneous Markov process. 
The diffusion X, starting at time s > 0, is given as the unique solution of 
the stochastic integral equation 

X t = X s + J [V(X U ) - b £ ' x °(u,X u )] du + y/i(W t - W s ), t > 8. 
By shifting the starting time back to the origin, this equation translates into 

X t+S = X S + f\v(X u+s ) - b £ > xo (u + s,X u+s )} du + y/iW t s , t > 0, 
Jo 

where W s is the Brownian motion given by Wf = Wt+ S — W s , which is 
independent of X s . Since we are mainly interested in the law of X, we may 
replace W s by W. 

For an initial condition £o £ an d s > 0, we denote by £ s, £° the unique 
solution of the equation 

(3.3) & = £o+ fvitJ-lfWfa + a^du + y/eWt, t>0. 

Jo 

Note that = X, and that £ s '£° has the same law as Xt+ S , given that 
X s = £,o- The interpretation of b e ' x ° as an interaction drift is lost in this 
equation, since b £,x ° does not depend on 

Now recall that b £ > x °(t, x) = E{<3?(x - Xf )}, which tends to <£(x - ^j(xo)) 
by Corollary 3.2. This motivates the definition of the following analogue of 
£ s '£°, in which u\ is replaced by the Dirac measure in ipt(xo). We denote by 
Y s ' y the solution of the equation 

(3.4) Y t = y+ [ V(Y U ) - ^(Y u - ^ +s (x )) du + y/iW t , t > 0. 

Jo 
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This equation is an SDE in the classical sense, and it admits a unique 
strong solution by Proposition 2.5. Furthermore, it is known that Y s ' y sat- 
isfies a large deviations principle in the space 

C 0T = {/ : [0, T] -> R d \f is continuous}, 

equipped with the topology of uniform convergence. This LDP describes 
the deviations of Y s,y from the deterministic system cpt = V(<pt) — &(<pt — 
ipt+s(xo)) with ipo = y. Observe that ip coincides with ip(xo) in case y = xq, 
and that nonexplosion of ip is ensured by the dissipativity properties of V 
and <E> as follows. By (2.4) we have 

— \\ipt - ipt+s\\ 2 =2{(p t - ipt+s,<Pt ~ ipt+s) 
(3.5) = 2{<pt - il>t+s,V&t) ~ ®{<Pt ~ i>t+s) ~ V(ip t+S )) 

<2(<pt-1>i+.,V(<Pt)-V(ik+s))- 

Since the last expression is negative for large values of \\ipt — ipt+s\\ by (2.8), 
this means that (pt — ipt+ s is bounded. But ip is bounded, so ip is also bounded. 

Let pvr(f,g) ■= sup < t < T ||/ - g\\ (f,ge C 0T ) be the metric correspond- 
ing to uniform topology, and denote by Hy the Cameron-Martin space of 
absolutely continuous functions starting at y that possess square-integrable 
derivatives. 



Proposition 3.3. The family (Y s>y ) satisfies a large deviations princi- 
ple with good rate function 

(3.6) I$(<p) = {i r 

oo, otherwise 

More precisely, for any closed set F C Cqt we have 



\ \\ip t -V(cpt) + $((pt-ipt+ s {xo))\\dt, if<p€H}j, 

J 



limsupelogP(y s ' s/ G F) < - inf 1$ (</>), 

and for any open set G C Cqt 

liminfelogP(Y^ G G) > - inf 1$ (</>). 

PROOF. Let a(t, y) := V(y) — $(y — tpt), and denote by F the function 
that maps a path g G Cqt to the solution / of the ODE 

ft = x + [ a(s,f s )ds + g t , 0<t<T. 
Jo 

Fix g G Cot, and let R > such that the deterministic trajectory ip(xo) as 
well as / = F(g) stay in -Br(O) up to time T. Note that nonexplosion of / is 
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guaranteed by dissipativity of a, much like in (3.5). Now observe that a is 
locally Lipschitz with constant 2K 2 r on -Br(O), uniformly w.r.t. t £ [0,T]. 
Thus, we have for g G Cqt-, f = F(g) such that / does not leave -Br(O) up 
to time T: 

\\ft - /t|| < 2K 2R f \\f s - f s \\ ds + \\g t -g t \\, 
Jo 

and Gronwall's lemma yields 

PoHf,f)<PoT(g,g)e 2K ™ T , 

that is, F is continuous. Indeed, the last inequality shows that we do not 
have to presume that / stays in -Br(O), but that this is granted whenever 
PoT(g,g) is sufficiently small. 

Since F is continuous and F(s/eW) = Y, we may invoke Schilder's theo- 
rem and the contraction principle, to deduce that Y satisfies a large devia- 
tions principle with rate function 

IotM = inf || £ \\gtfdt : g G H^,F(g) = p j. 
This proves the LDP for (Y s ^). □ 

Notice that the rate function of Y measures distances from the deter- 
ministic solution ip just as in the classical case without interaction, but the 
distance of from ip is weighted by the interaction between the two paths. 

By means of the rate function, one can associate to Y s ' v two functions 
that determine the cost, respectively energy, of moving between points in 
the geometric landscape induced by the vector field V. For t > 0, the cost 
function 

C s (y,z,t)= inf 4f(f), y,zeR d , 
fec f.ft=x 

determines the asymptotic cost for the diffusion Y s ' y to move from y to z 
in time t, and the quasi-potential 

Q s (y,z)=mW s (y,z,t) 
describes its cost of going from y to z eventually. 

3.2. Large deviations principle for the self-stabilizing diffusion. We are 
now in a position to prove large deviations principles for £ and X by showing 
that £ and Y are close in the sense of large deviations. 
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Theorem 3.4. For any e > let Xg,£o £ IR d converge to some xq G 
R d , respectively y G M d , as e — > 0. Denote by X £ the solution of (3.1) starting 
at Xq. Let s > 0, and denote by ft i/ie solution of (3.3) starting in ft mi/i 
time parameter s, that is, 

(3-7) ft = ft + f + + M t>0, 

waere 6 e ' x °(t,x) =E[$(x - Xf)]. 

Taen i/ie diffusions (ft) e >o satisfy on any time interval [0, T] a Zarge 
deviations principle with good rate function (3.6). 

Proof. We shall show that £ := ft is exponentially equivalent to Y := 
Y s ' v as defined by (3.4), which has the desired rate function; that is, we 
prove that for any 5 > we have 

(3.8) limsupelogP(por(£,JO >S) = -oo. 

e^O 

Without loss of generality, we may choose R > such that Xq, y G Br(0) and 
that tpt(xo) does not leave -Br(O) up to time s + T, and denote by er# the 
first time at which £ or Y exits from Br(0). Then for t < or 

lift - it|| < lift - 2/11 + f l|v(&) - K(y„)|| du 

(3.9) ( i0 

+ / \\b £ > x o(u + s,£ u )-$(Y u -i> u+s (x ))\\du. 





The first integral satisfies 

f 11^(60 - v(y u )\\ du < k r f ||ft - r u || d«,t < *r, 

Jo Jo 

due to the local Lipschitz assumption. Let us decompose the second integral. 
We have 

\\b £ ' x °(u + s,ft) - <S>(Y U - ^ u+s (x ))|| 

<||6^o( u + S) ft)-$(ft-^ +s (4))|| 

+ ||$(ft - ^u+s(x £ )) ~ $(ft - 1pu+s(x ))\\ 

+ ||$(ft - ijj u+s (x )) - <S>(Y U - i; u+s (x ))\\. 

Bounds for the second and third term in this decomposition are easily de- 
rived. The last one is seen to be bounded by -f^flllft — Yu\\> since ft Y as 
well as ip are in Br(0) before time gr A T. For the second term we also use 
the Lipschitz condition to deduce that 

||$(ft - ip u+s (x e Q )) - $(ft - ip u+s (x ))\\ < K 2 r\\^ u + s {xI) - ip u+s (x )\\. 
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As a consequence of the flow property for tp this bound approaches as 
e — > uniformly w.r.t. u E [0,T]. 

By combining these bounds and applying Gronwall's lemma, we find that 

Ut-Y t \\ 

(3.10) <exp{2K2 R t}(\\So-y\\ + K 2R J \\ip u+s (x £ )-ip u+s {x )\\du 

+ J \\b £ ' x Hu + s,£ u ) - *(f u - ^u+s (x £ ) ) 1 1 duj 

for t < (Jr. Since £ is bounded before an the r.h.s. of this inequality tends 
to zero by Corollary 3.2. 

The exponential equivalence follows from the LDP for Y as follows. Fix 
5 > 0, and choose £o > such that the r.h.s. of (3.10) is smaller than 6 for 
£ < £o- Then \\^ t — Y t \\ > 6 implies that at least one of £t or Y t is not in B R (0), 
and if ^ Br(0), then Yi ^ fi^/ 2 (0) if 5 is small enough. Thus we can bound 
the distance of £ and Y by an exit probability of 1". For I > let 77 denote 
the diffusion Y's time of first exit from 5/(0). Then, by Proposition 3.3, 

limsupelogP(p T(£, Y) > 5) 

(3.11) < limsupelogP(r R/2 < T) 

e^O 

< -ini^C s (y,z,t):\z\ > -,0<t<rj. 
The latter expression approaches — 00 as R — > 00. □ 

Theorem 3.4 allows us to deduce two important corollaries. A particular 
choice of parameters yields an LDP for X, and the e-dependence of the 
initial conditions permits us to conclude that the LDP holds uniformly on 
compact subsets, a fact that is crucial for the proof of an exit law in the 
following section. The arguments can be found in [7]. 

Let F Xo (X £ ") denote the law of the diffusion X starting at xq S ffi d . 

Corollary 3.5. Let L C M, d be a compact set. For any closed set F C 
Cot we have 

limsupelog sup F Xo (X &F)<- inf inf 

e~>0 x £L x £L<j>&F 

and for any open set G C Cqt 

liminfelog inf F Xo (X e G) > - sup inf ^W- 

x GL XoG L<t>£G 
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Proof. Choosing Xq = £q an d s = implies £ e = X s in Theorem 3.4, 
which shows that X satisfies an LDP with rate function Iqt ° ■ Furthermore, 
this LDP allows for e-dependent initial conditions. This implies the uni- 
formity of the LDP, as pointed out in the proofs of Theorem 5.6.12 and 
Corollary 5.6.15 in [7]. Indeed, the e-dependence yields for all xq E M. d 

limsup elogF y (X E F) < - inf (0), 

for otherwise one could find sequences e n > and y n E R rf such that e n — ► 0, 
y n — > x and 

limsupe^logPjJ* E F) > - inf (</>). 

n— >oo <j>EF 

But this contradicts the LDP. 

Now the uniformity of the upper large deviations bound follows exactly 
as demonstrated in the proof of Corollary 5.6.15 in [7]. The lower bound is 
treated similarly. □ 



The next corollary is just a consequence of the e-dependent initial condi- 
tions in the LDP for £. 



Corollary 3.6. Let L C M. d be a compact set. For any closed set F C 
Cot we have 

limsupelog sup F(£ s ' x ° E F) < - inf inf (0), 

e^O x eL x £L<f>£F 

and for any open set G C Cot 

liminf elog inf P(£''*° E G) > - sup inf Iff {<)>). 

£^0 x £L XQ£ L<i>£G 



3.3. Exponential approximations under stability assumptions. The aim 
of this subsection is to exploit the fact that the inhomogeneity of the diffu- 
sion Y s ' y is weak in the sense that its drift depends on time only through 
ipt+s(xo). If the dynamical system ip = V(i/j) admits an asymptotically sta- 
ble fixed point x s t a bie that attracts xq, then the drift of Y s ' y becomes almost 
autonomous for large times, which in turn may be used to estimate large 
deviations probabilities for ^ s,y . We make the following assumption. It will 
also be in force in Section 4, where it will keep us from formulating results 
on exits from domains with boundaries containing critical points of DV, in 
particular saddle points in the potential case. 
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Assumption 3.7. 

(i) Stability: there exists a stable equilibrium point x s tabic G ^ d of the 
dynamical system 

(ii) Convexity: the geometry induced by the vector field V is convex, 
that is, the condition (2.6) for V holds globally: 

(3.12) (h,DV(x)h) <-K v 
for h£R d s.t. \\h\\ = 1 and x G R d . 

Under this assumption it is natural to consider the limiting time-homogeneous 
diffusion Y°°' y defined by 

(3.13) dY t °° = V(Y t °°) dt - $(Y t °° - x stable ) dt + ^e dW t , Y °° = y. 

Lemma 3.8. Let L C M. d be compact, and assume that actable attracts all 
y G L, that is, 

lim tp t (y) = Stable Vy G L. 

t— >oo 

Then Y°°' y is an exponentially good approximation ofY s,y , that is, for any 
5 > we have 

lim lim sup e log sup P(p 0T (Y s ' y , Y°°' y ) > 5) = -oo. 

r ^°° e^O y£L,s>r 

Proof. We have 

llY s,y _ y ^ < I* m ys, y) _ V( Y^V)\\ du 
JO 

+ f\MY^ y - Vw(y)) - HY™« - XstMe )\\ du. 

Jo 

Let a S jf be the first time at which Y s ' y or Y°°' y exits from Br(0). For 
t < ^jfi we may use the Lipschitz property of 3> and V, to find a constant 
cr > s.t. 

||y.,« _ yoo|| < CR f* ||y., tf _ yoo|| ^ + Cfl Tp T (^+.(2/), Stable)- 
JO 

By assumption the second term converges to as s — > oo, uniformly with 
respect to y G L since the flow is continuous with respect to the initial data. 
Hence, by Gronwall's lemma there exists some r = r(R,5) > such that for 

s > r 

sup sup \\Y t s ' y -Y t °°\\ < 5/2. 

y<ELo<t<a a ' y 
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We deduce that 



¥(p 0T (Y s ' y , y°°) > 5/2) < V(t» < T) Vs > r, y € L, 



where for I > rf denotes the first exit time of Y°°' y from B[(0). Sending 
r, R — > oo and appealing to the uniform LDP for Y 00 ^ finishes the proof, 
much as the proof of Theorem 3.4. □ 

This exponential closeness of Y°°' y and Y s ' y carries over to ^ s ' y under 
the aforementioned stability and convexity assumption, which enables us 
to sharpen the exponential equivalence proved in Theorem 3.4. In order to 
establish this improvement, we need a preparatory lemma that strengthens 
Corollary 3.2 to uniform convergence over the whole time axis. This unifor- 
mity is of crucial importance for the proof of an exit law in the next section 
and depends substantially on the strong convexity assumption (3.12). 

Lemma 3.9. We have 



Proof. Let f(t) := E(||Z t || 2 ), where Z t = X t - ^(x ). In the proof of 
Lemma 3.1 we have seen that 

f'{t) < 2E[(Z t , V(Z t + VtOr )) " V(Mxo)))} + ed 
< -2K v E(\\Z t \\ 2 ) + ed = -2K v f(t) + ed. 

This means that {t > 0:f'{t) < 0} D {t > 0:/(t) > ^7}. Recalling that 

/(0) = 0, this allows us to conclude that / is bounded by 27^7- Now an 
appeal to the proof of Corollary 3.2 finishes the argument. □ 

Proposition 3.10. Let LcR d be compact, and assume that x s t a bie 
attracts all y £ L. Then Y°°' y is an exponentially good approximation of 
£ s > y , that is, for any 5 > we have 



Proof. Recall the proof of Theorem 3.4. For y E L and s > we have 




lim limsupelog sup F(p 0T (£ s > y , Y°°' y ) > 5) 



e— >0 y&L,s>r 



= — OO. 



(3.14) 
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for t < which denotes the first time that £* ,v or Y s,y exits from Br(0). 
By Lemma 3.9, the integrand on the r.h.s. converges to zero as e — > 0, uni- 
formly w.r.t. s > 0. Therefore, if we fix 5 > 0, we may choose R = R(S) 
sufficiently large and Eq > such that for e < Eq, and all s > 

P(jHxr(t'> y ,Y'»)>$)<F(T% 2 <T) 

< P(r^'| < T) + P(p r(^ oo ' y , F^) > 

where for I > 0, < s < oo, tv denotes the first exit time of the diffusion 
y s 'f from the ball 5/(0). By the uniform LDP for Y°°' v and Lemma 3.8 the 
assertion follows. □ 

4. The exit problem. As a consequence of the large deviations principle, 
the trajectories of the self-stabilizing diffusion are attracted to the determin- 
istic dynamical system tp = V(ip) as noise tends to 0. The probabilities of 
deviating from tp are exponentially small in e, and the diffusion will certainly 
exit from a domain within a certain time interval if the deterministic path 

exits. The problem of diffusion exit involves an analysis for the rare event 
that the diffusion leaves the domain although the deterministic path stays 
inside, that is, it is concerned with an exit which is triggered by noise only. 
Clearly, the time of such an exit should increase as the noise intensity tends 
to zero. In this section we shall derive the precise large deviations asymp- 
totics of such exit times, that is, we shall give an analogue of the well-known 
Kramers-Eyring law for time- homogeneous diffusions. 

Let us briefly recall this law, a detailed presentation of which may be found 
in Section 5.7 of [7]. For further classical results about the exit problem we 
refer to [5, 6, 8] and [17]. 

A Brownian particle of intensity e that wanders in a geometric landscape 
given by a potential U is mathematically described by the classical time- 
homogeneous SDE 

dZf = -VU{Zf)dt + VsdWt, Z^ = x €R d . 

If x* is a stable fixed point of the system x = — W{x) that attracts the initial 
condition xq and t £ denotes the exit time from the domain of attraction of 
x* , then the asymptotics of t 6 is described by the following two relations: 

(4.1) limelogE(r £ ) =U, 

(4.2) lim F(e {0 ~ 5 ^ £ <t £ < e i0+5)/e ) = 1 V£ > 0. 

Here U denotes the energy required to exit from the domain of attraction 
of x* . This law may roughly be paraphrased by saying that t £ behaves like 
exp y as e — > 0. 



KRAMERS' TYPE LAW 



33 



Let us now return to the self-stabilizing diffusion X s , defined by (3.1). 
Intuitively, exit times should increase compared to the classical case due to 
self-stabilization and the inertia it entails. We shall show that this is indeed 
the case, and prove a synonym of (4.1) and (4.2) for the self-stabilizing 
diffusion. Our approach follows the presentation in [7]. 

Let D be an open bounded domain in M. d in which X s starts, that is, 
xq £ D, and denote by 

t £ d = mf{t > : X s £ dD} 

the first exit time from D. We make the following stability assumptions 
about D. 



Assumption 4.1. 

(i) The unique equilibrium point in D of the dynamical system 

(4.3) ^ = Vtyt) 

is stable and given by x sta bic € D. As before, ipt{xo) denotes the solution 
starting at xq. We assume that lim^oo ipt{%o) = ^stable- 

(ii) The solutions of 

(4.4) (f> t = V((pt) ~ $(<f>t - ^stable) 

satisfy 

4> eD <f>t£D Vi>0, 

4>0 e D =>• lim (j) t = X s t a ble- 

The description of the exponential rate for the exit time of ltd diffusions 
with homogeneous coefficients was first proved by Freidlin and Wentzell via 
an exploitation of the strong Markov property. The self-stabilizing diffusion 
X s is also Markovian, but it is inhomogeneous, which makes a direct appli- 
cation of the Markov property difficult. However, the inhomogeneity is weak 
under the stability Assumption 4.1. It implies that the law of Xf converges 
as time tends to infinity, and large deviations probabilities for X £ may be 
approximated by those of Y°° in the sense of Proposition 3.10. Since Y°° 
is defined in terms of an autonomous SDE, its exit behavior is accessible 
through classical results. The rate function that describes the LDP for Y°° 
is given by 

(4.5) I$ v ((p) = I \J Q ~ V ^ + $ (^* _a; stabie) II 2 dt, i£<peH$, 

[ oo, otherwise. 

The corresponding cost function and quasi-potential are defined in an ob- 
vious way and denoted by C°° and respectively. The minimal energy 
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required to connect the stable equilibrium point arable to the boundary of 
the domain is assumed to be finite, that is, 

Qoo := inf Q°°(> stable ,z) < oo. 
z£dD 

The following two theorems state our main result about the exponential 
rate of the exit time and the exit location. 

Theorem 4.2. For all xo £ D and all rj > 0, we have 

(4.6) limsupelog{l -P xo (e ( ^-- r?)/£ <t £ d < e (Q- +r ?)/ £ )} < - v /2 

and 

(4.7) limelogE a!0 (7^) = Q oo . 

e— >0 

Theorem 4.3. If N C dD is a closed set satisfying 

inf r Q°°{x stah ie,z) > Qoo> 

then it does not see the exit point: for any xq G D 

limP x , ) (^ D GiV) = o. 

The rest of this section is devoted to the proof of these two theorems. 
In the subsequent section, these results are illustrated by examples which 
show that the attraction part of the drift term in a diffusion may completely 
change the behavior of the paths, that is, the self-stabilizing diffusion stays 
in the domain for a longer time than the classical one, and it typically exits 
at a different place. 

4.1. Enlargement of the domain. The self-stabilizing diffusion lives in 
the open, bounded domain D which is assumed to fulfill the previously 
stated stability conditions. In order to derive upper and lower bounds of 
exit probabilities, we need to construct an enlargement of D that still enjoys 
the stability properties of Assumption 4.1(h). This is possible because the 
family of solutions to the dynamical system (4.4) defines a continuous flow. 

For 5 > we denote by D s := {y G M. d : dist(y, D) < 6} the open <5-neighbor- 
hood of D. The flow (j> is continuous, hence uniformly continuous on D due 
to boundedness of D, and since the vector field is locally Lipschitz. Hence, if 
5 is small enough, the trajectories 4>t(y) converge to actable f° r V € D s , that 
is, for each neighborhood V C D of x s t a bie there exists some T > such that 
for y G D s we have (f>t(y) G V for all t>T. Moreover, the joint continuity of 
the flow implies that, if we fix c > 0, we may choose 5 = 5(c) > such that 

sup{dist (0 t (y) , D) : t G [0, T] , y G D S } < c. 
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Let 

O s = {y£R d : sup dist(0t(y),-D) < c,(fyp(y) € v\. 
L te[o,T] ) 

Then O s is a bounded open set which contains D s and satisfies Assumption 
4.1(h). Indeed, if 8 is small enough, the boundary of O s is not a characteristic 
boundary, and n<5>o O s = D. 

4.2. Proof of the upper bound for the exit time. For the proof of the two 
main results, we successively proceed in several steps and establish a series of 
preparatory estimates that shall be combined afterward. In this subsection, 
we concentrate on the upper bound for the exit time from D, and establish 
inequalities for the probability of exceeding this bound and for the mean 
exit time. 

In the sequel, we denote by F StV the law of the diffusion ^ s,y , defined 
by (3.3). Recall that by the results of the previous section, £ s ' y satisfies a 
large deviations principle with rate function I s ' y . The following continuity 
property of the associated cost function is the analogue of Lemma 5.7.8 in 
[7] for this inhomogeneous diffusion. The proof is omitted. 

Lemma 4.4. For any 5 > and s £ [0,oo), there exists g > such that 

(4.8) sup inf C s {x,y,t)<5 

x,y&B e (x stMc )t£[0M 

and 

(4.9) sup inf C s (x,y,t) < 5, 

(x,2/)er*e[o,i] 

where F = {(x, y) : mf zedD (\\y - z\\ + \\x - z\\) < g}. 

Let us now present two preliminary lemmas on exit times of £ s ' y . In slight 
abuse of notation, we denote exit times of £ s,y also by rfy, which could 
formally be justified by assuming to look solely at the coordinate process 
on path space and switching between measures instead of processes. On 
the other hand, this notation is convenient when having in mind that £ s > y 
describes the law of X s restarted at time s, and that X s may be recovered 
from £ s ' y for certain parameters. 

Lemma 4.5. For any r\ > and g > small enough, there exist Tq > 0, 
so > and £q > such that 

for all e < £o and s > s$. 

y&B e (x BtMc ) 
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Proof. Let g be given according to Lemma 4.4. The corresponding re- 
sult for the time-homogeneous diffusion Y°° ,y is well known (see [7], Lemma 
5.7.18), and will be carried over to £ s ' y using the exponential approximation 
of Proposition 3.10. Let P^y denote the law of Y°°' y . The drift of Y°°> y is lo- 
cally Lipschitz by the assumptions on V and and we may assume w.l.o.g. 
that it is even globally Lipschitz. Otherwise we change the drift outside a 
large domain containing D. 

If 5 > is small enough such that the enlarged domain O s satisfies As- 
sumption 4.1(h), Lemma 5.7.18 in [7] implies the existence of £\ and To such 
that 

(4.10) inf Poo,j/(^ < T ) > e-^-^/ 3 ^ foralle<£i. 

y£B e (x sta , hlc ) 

Here denotes the minimal energy 

QL = inf E Q°° (^stable, z). 
z£dO s 

The continuity of the cost function carries over to the quasi-potential, that 

is, there exists some So > such that IQoo — Qoo\ < ??/3 for 5 < <5o- 

Now let us link the exit probabilities of Y°°' y and ^ s ' y . We have for s > 

K,v(r E D < T ) 

(4.11) > P({^ y exits from D before T } n {po,T (t s ' y ,Y°°' y ) < 5}) 
>Poc,y(r^ <T ) -P^To^,^ 00 ) > 6). 

Moreover, by the exponential approximation we may find £2 > and sq > 
such that 

sup P(p 0j T (£ s ' y \Y°°) >5)< e -^ +v/2)/£ VS>S ,£<£2- 
y&B g (x Bt!lblc ) 

Since D s C O s , we deduce that for e < £q = E\ A £2 and s > sq 

inf P s , y ( T b < T o) > e~^ +ri/3)/£ - e -^ + ^ /2)/e > e -^ +v)/£ . 

y&B e (x stMo ) ' LI 

By similar arguments, we prove the exponential smallness of the prob- 
ability of too long exit times. Let S e = inf{i > 0:^'^ £ -B e (x sta bic) U dD}, 
where g is small enough such that B g (x s tahie) is contained in the domain D. 

Lemma 4.6. For any g > sufficiently small and for any K > there 
exist eq > 0, T\ > and r > such that 

sup P St y(T, e > t) < e- K/e Vt>Ti. 

y&D,s>r 
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Proof. As before, we use the fact that a similar result is already known 
for Y°°' y . For 5 > small enough, let 

= M{t > : Y t °° 6 B^Cxstabie) U aO 5 }. 

By Lemma 5.7.19 in [7], there exist T\ > and ei > such that 

su P P 00i j / (S^>t)<e-^ /£ Vt>Ti;e<ei. 



Now the assertion follows from 



supP s ,j,(E e > Ti) < supPoo^CEj > Ti) + supP(po,T 1 (^ ,y ,i" 30,y ) > S), 

yeD yeD yeD 

since the last term is exponentially negligible by Proposition 3.10. □ 

The previous two lemmas contain the essential large deviations bounds 
required for the proof of the following upper bound for the exit time of X £ . 

Proposition 4.7. For all x e D and rj > we have 

(4.12) limsupe logP^rf, > e^^ +v)/£ ) < -i]/2 

£->0 

and 

(4.13) lim sup e log E Xo [t £ d ] < . 

Proof. The proof consists of a careful modification of the arguments 
used in Theorem 5.7.11 in [7]. By Lemmas 4.5 and 4.6, there exist T = 
To + T\ > 0, £o > and ro > such that for T > T, e < Eq and r > tq we 
have 

^:=infP r , tf (r&<r) 

yeD 

(4.14) > inf P r|2/ (S, < Ti) inf P s>y (r|, < T ) 

2/€-° 2/€B e (a; atab i e ),s>r 

> exp {-^±^} =; ,». 

Moreover, by the Markov property of £ s ' v , we see that for fcsN 
P S0 (r& > 2(k + 1)T) = [1 - F xo (t e d < 2(k + l)T\r £ D > 2kT)\ 

xP XQ (T £ D >2kT) 

< 1 - inf F 2kT , y (t £ d < 2T)\ P X0 (t £ d > 2kT) 
yen j 

<(l- q 2 2 !f)F X0 (r £ D >2kT), 
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which by induction yields 



fc-i 



(4.15) 



■i-t) 



(Th>2kT)<Yl(l-q 2 2 F). 



i=Q 



Let us estimate each term of the product separately. We have 
l-g£r =sup¥ 2i T,y(r!j>2T) 

y&D 

< sup Par^Cr!, > T) supP (2i+1)T . ?/ (r|, > T) 

yGD y&D 

<sup¥ {2l+my (r £ D >T). 

y&D 

By choosing T large enough, we may replace the product in (4.15) by a 
power. Indeed, for T > max(T, tq) we have (2i + 1)T > ro for all i G N, 
which by (4.14) results in the uniform upper bound 

i-<z!f <i-# +1)T <i-^. 

By plugging this into (4.15), we obtain a "geometric" upper bound for the 
expected exit time, namely 



K X0 [t £ d }<2T 



l + ^supP X0 (r| ) >2A : r) 



< 2T 



< 2T 



oo k—1 
fc=l i=0 



1 + E( 1 -9t)" 



fe=i 



2T 
9^" 



This proves the claimed asymptotics of the expected exit time. Furthermore, 
an application of Chebyshev's inequality shows that 



c (T e D >e {Q ~> +r > )/£ )< 



e x [t £ d ] <2T £. 



qf 



which is the asserted upper bound of the exit probability. □ 



2Te- v/2e , 



4.3. Proof of the lower bound for the exit time. In order to establish the 
lower bound of the exit time, we prove a preliminary lemma which estimates 
the probability to exit from the domain D \ .B e (actable) at the boundary of 
D. This probability is seen to be exponentially small since the diffusion is 
attracted to the stable equilibrium point. Let us denote by S g the boundary 
of B g (x sta bi e ), and recall the definition of the stopping time T, g . 
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Lemma 4.8. For any closed set N C dD and r\ > 0, there exist Eq > 0, 
go > and ro > such that 

elog sup P(gj» €N)<- inf Q°°(x stablc , z)+n 
for all e < Eq, r>r and g < g . 

Proof. For 5 > we define a subset S s of D s by setting 
S 5 ■= D s \{yeR d : dist(y, N) < 5}. 

Furthermore, let 

M 5 := OS 5 n{y€R d : dist(y, N) < 5}. 

S s contains the stable equilibrium point arable; an d as such it is unique in 
S s if 5 is small enough. 



The proofs of Lemma 5.7.19 and Lemma 5.7.23 in [7] can be adapted to 
the domain S s , since an exit of the limiting diffusion Y°° from the domain 
O s defined in Section 4.1 always requires an exit from S s . Hence, there exist 
Ei > and gi > such that 

elog sup F^yiY^eM 5 ) < - inf .Q°°0r stablc ,z) + \ 

for e < Ei and g < gi, where denotes the first exit time from the domain 
S s \ ^(^stabic)- By the continuity of the quasi-potential, we have 

T] 

- inf Q°°(a; sta bie,2) + - < ~ inf Q°° (^stable, z) +V 

if 6 > is small enough. Therefore, it is sufficient to link the result about 
the limiting diffusion to the corresponding statement dealing with £ S,J/ . By 
Lemma 4.6, we can find Ti > 0, e\ > and r\ > such that 

elog sup P S)2/ (£ e > Ti) 

y£S2 g ,s>r 

(4-16) 

< - inf Q°° (a; sta bie, z) + 77 Ve<ei,r>ri. 
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If £ e < T x and p , Tl (C' y , Y°°) < 5, then {gj G N} is contained in {Y^ G 
A^}. Thus, 

G iV) < P(gJ G N, E c < Tx) + P Si3/ (S e > Tx) 

< p(y^ g aa- 5 ) + p(po >Tl Y°°>y) > 5) 

Q 

+ p Sj?/ (s e >r 1 ). 

By (4.16) and Proposition 3.10, the logarithmic asymptotics of the sum on 
the r.h.s. is dominated by the first term, that is, the lemma is established. 

□ 

We are now in a position to establish the lower bound for the exit time 
which complements Proposition 4.7 and completes the proof of Theorem 4.2. 

Proposition 4.9. There exists rjo > such that for any r\ < rjo 

(4.17) limsupelogP^frf) < e ( ^°°^ )/e ] < -ry/2 

and 

(4.18) limmfelogE^rf)]^. 

Proof. In a first step we apply Lemma 4.8 and an adaptation of 
Lemma 5.7.23 in [7]. The latter explains that the behavior of an ltd dif- 
fusion on small time intervals is similar to the behavior of the martingale 
part, which in our situation is simply given by y/eWt- We find ro > 0, T > 
and £o > such that for e < e$ 

sup P(gJ G 3D) < e -(Qoo-W2)/ £) 

y£S2 e ,s>ro 

sup f( sup \\& y - y\\ >q]< e-^^ /2)/e . 

y£D,s>r \0<t<T ) 

In the sequel, we shall proceed as follows. First, we wait for a large 
period of time r\ until the diffusion becomes "sufficiently homogeneous," 
which is possible thanks to the stability assumption. Since x s tabic attracts 
all solutions of the deterministic system, we may find r\ > ro such that 
Vv(^o) G -Bg(x s tabie) fc> r r > r i- Second, after time r±, we employ the usual 
arguments for homogeneous diffusions. Following [7], we recursively define 
two sequences of stopping times that shall serve to track the diffusion's ex- 
cursions between the small ball -B e (actable) around the equilibrium point and 



KRAMERS' TYPE LAW 



41 




the larger sphere = 9B2 g (actable); before it finally exits from the domain 
D. 

Set $o = f\ , and for m > let 

r m = inf {t > m : 6B e U <9L>} 

and 

$ m+1 = inf{t > r m : Xf G S^}. 
Let us decompose the event {rfj < kT + r{\. We have 

Vx (iD < kT + ri) < W X0 {{t £ d < n} U {X* £ J3 2B (a; B tabie)}) 

(4.19) 

+ sup P a , y (ijj < kT). 

y&S 2g ,s>r 1 

The first probability on the r.h.s. of this inequality tends to as e — > 0. 
Indeed, by the large deviations principle for X s on the time interval [0, n], 
there exist % > and e 2 > such that 

elog¥ Xo ({r £ D < n} U £ £ 2£ >(x stablc )}) < -r?/2 

for e < £2 and r/ < 7/0 ■ For the second term in (4.19), we can observe two 
different cases: either the diffusion exits from D during the first k exits from 
D \ Bg(x 8 tab\e), ° r the minimal time spent between two consecutive exits is 
smaller than T. This reasoning leads to the bound 

W s ,y(Th <kT)<J2 Vs,y( T b = T m ) + f>s,y ( min (0 TO - T m _i) < T 

\l<m<k 

m=0 \ — — 

Let us now link these events to the probabilities presented at the beginning 
of the proof. We have 

sup P S)2/ (rB = r m )< sup F 8tV (£ v edD) 
yeS2 B ,s>ri yeS2 e ,s>ro 
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and 

sup P S)3/ ((i? m - r m -i) < T) < sup W S J sup \\£t V -y\\>Q 
y&S 2e ,s>n yeS2 B ,s>r \0<t<T 

which yields the bound 

sup P, )W (7£ < kT) < (2k + lje-^--"/ 2 )/ 6 . 

y£S2 e ,s>ri 

Thus, by choosing A; = [(e^°° _,? ^ e — n)/Tj + 1, we obtain from (4.19) 

that is, (4.17) holds. Moreover, by using Chebyshev's inequality, we obtain 
the claimed lower bound for the expected exit time. Indeed, we have 

® X0 (T e D ) > e^~-")/ £ (l - P xo (rh < e®~~^)) 

> e (5 00 -'»)A( 1 _ (i + ST^e"^), 

which establishes (4.18). □ 

We end this section with the proof of Theorem 4.3 about the exit location. 

Proof of Theorem 4.3. We use arguments similar to the ones of the 
preceding proof. Let 

Q oc (N) = inf Q°°(x stable ,z), 

and assume w.l.o.g. that Q 00 (N) < oo. Otherwise, we may replace Qoo(N) 
in the following by some constant larger than . As in the preceding proof, 
we may choose T > 0, tq > and Eq > such that 

sup P., v (gi v G cW) < e-^-W-^A Ve < £o> 
sup P s>2/ f sup ||& s,2/ - >q)< e-&^- T il 2 V £ Ve < e . 

y£D,s>r Q ' \0<t<T J 

It suffices to study the event A = {rfj < kT + ro} n G N} for positive 
integers fc. We see that 

P x . (A) < P X0 (X £ Q $ B 2e (x staUc )) + sup P Sty {T £ D < kT) 

y£S2 e ,s>ro 
k 

< P Xo (X £ i B 2e (x stahlc )) + Vs,y(rh = T m ,C'i G N) 

+ P S ,J min (i9 m - r m _i) < T 

< P X0 (A^ i B 2g (x staUe )) + (2k + l) e -<5~ W-^/2)A. 
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The choice k = [(e^o.W-v)^ _ rQ )/T\ + 1 yields 

F X0 (A) < F xo (X' o <£ B 2Q (x stahlc )) + 5T- 1 e-"/ 2£ . 

This implies that F Xo {t £ d < e^W-f)/^^ £ iV) — > as e — > 0. Now choose 77 
small enough such that Q oc (A r ) — 77 > + r/. Then Proposition 4.7 states 
that the exit time of the domain D is smaller than e^°° +ri ^' £ with probability 
close to 1. The combination of these two results implies ¥ Xo (X £ ^ £ N) — > 
as e — » 0. □ 

5. The gradient case: examples. The structural assumption about 
namely its rotational invariance as stated in (2.4), implies that <J> is always 
a potential gradient. In fact, this assumption means that $ is the gradient 
of the positive potential 

r\\x\\ 

A(x) = / 4>(u) du. 
Jo 

In this section, we make the additional assumption that the second drift com- 
ponent given by the vector field V is also a potential gradient, which brings 
us back to the very classical situation of gradient type time-homogeneous 
Ito diffusions. In this situation, quasi-potentials and exponential exit rates 
may be computed rather explicitly and allow for a good illustration of the 
effect of self-stabilization on the asymptotics of exit times. 

We assume from now on that V = — VJ7 is the gradient of a potential U 
on R d . Then the drift of the limiting diffusion Y°° defined by (3.13) is also 
a potential gradient, that is, 

b(x) := V(x) - $(x - Stable) = - V(C7 (x) + A(x - Stable))- 

A simple consequence of Theorem 3.1 in [8] allows one to compute the quasi- 
potential explicitly in this setting. 

Lemma 5.1. Assume that V = —VU. Then for any z € D, 

Q°° (^stable, z) = 2(U(z) -U(x sta ble) + A(z - Stable))- 

In particular, 

Qoo = inf 2(U(z) -U{x stah i e ) + A(z - arable))- 

Observe that the exit time for the self-stabilizing diffusion is strictly larger 
than that of the classical diffusion defined by 



dZf = V{Z e t )dt + ^edW t , Zl = x Q . 
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Indeed, by the theory of Preidlin and Wentzell, 

(5.1) lim e log E X0 (r £ D (Z £ )) = inf 2(U(z) - U(x 8taUe )) 

e— >0 z£oD 

(5.2) <Q oo = limelogE :C0 (r| ) (X e )). 

e— »0 

The exit problem is in fact completely different if we compare the dif- 
fusions with and without self-attraction. We have already seen that the 
exponential rate is larger in the attraction case. Let us next see by some 
examples that the exit location may change due to self-stabilization. 



5.1. The general one- dimensional case. In this subsection we confine 
ourselves to one-dimensional self-stabilizing diffusions. In dimension one, 
the structural assumptions concerning $ and V are always granted, and we 
may study the influence of self-stabilization on exit laws in a general setting. 
Let a < < b, and assume for simplicity that the unique stable equilibrium 
point is the origin 0. Denote by U(x) = — Jq V(u) du the potential that in- 
duces the drift V. As seen before, the interaction drift is the gradient of 
the potential A(x) = 4>{u) du. Since we are in the gradient situation, 
the exponential rate for the mean exit time from the interval [a, b] can be 
computed explicitly. 

If we denote by t x (X £ ) = inf{t > : Xf = x} the first passage time of the 
level x for the process X s and tj = r a A r^, then the exit law of the classical 
diffusion Z £ (i.e., without self-stabilization) is described by 

limP (e( Q "-^ < r 7 (Z e ) < e ^^l £ ) = 1 

and 

limelogEo(r / (Z e )) = Q-, 

e^O 

where Qo° = 2min([/(a), U{b)). Moreover, if we assume that U(a) < U(b), 
then P (t/(Z £ ) = T a {Z £ )) -> 1 as e -> 0. 

The picture changes completely if we introduce self-stabilization. The 
quasi-potential becomes 

Qf 5 = 2min([/(a) + A(a), U(b) + A(b)) > 

so the mean exit time of X s from the interval / is strictly larger compared 
to that of Z £ . This result corresponds to what intuition suggests: the process 
needs more work and consequently more time to exit from a domain if it 
is attracted by some law concentrated around the stable equilibrium point. 
Furthermore, if a and b satisfy 



A{b)-A{a)<U(a)-U(b), 
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we observe that ¥q(ti(X £ ) = Tb(X £ )) — > 1, that is, the diffusion exits the 
interval at the point b. Thus, we observe the somehow surprising behavior 
that self-stabilization changes the exit location from the left to the right 
endpoint of the interval. See Figure 1. 

5.2. An example in the plane. In this subsection, we give another explicit 
example in dimension two, in order to illustrate changes of exit locations in 
more detail. 

Let V = —VU, where 

U(x,y) = 6x 2 + y 2 , 
and let us examine the exit problem for the elliptic domain 
D = {(x,y)£R 2 :x 2 + ly 2 <l}. 

The unique stable equilibrium point is the origin actable = 0. 

The asymptotic mean exit time of the diffusion Zf starting in is given 
by lim e ^oelogEo(' r D(-^ £ )) = srn ce the minimum of the potential on dD 
is reached if y = ±2 and x = 0. Let us now focus on its exit location, and 
denote N^ x ^ = dD f] B g ((x,y)). The diffusion exits asymptotically in the 
neighborhood -/V"(q 2 ) with probability close to 1/2 and in the neighborhood 
^V(o,-2) with the same probability. 

Now we look how self-stabilization changes the picture. For the interaction 
drift we choose Q(x,y) = VA(x,y), with A(x,y) = 2x 2 + 2y 2 . First, the self- 
stabilizing diffusion X s starting in needs more time to exit from D, namely 
lim £ ^o £ l°gEo(r| ) (X £ )) = 16. More surprisingly, though, the exit location 
is completely different. The diffusion exits asymptotically with probability 
close to 1/2 in the neighborhoods iV(_ lj0 ) and -/V(i,o)> respectively. 




Fig. 1. Potentials U (left picture) and U + A (right picture). 
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Exit locution of the 
tJussical diffusion 



&e.l f-staMlizing diffusion 
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